Many-body effect in spin dephasing in n-typed GaAs quantum wells 
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By constructing and numerically solving the kinetic Bloch equations we perform a many-body 
study of the spin dephasing due to the D'yakonov-Perel' effect in n-typed GaAs (100) quantum wells 
for high temperatures. In our study, we include the spin-conserving scattering such as the electron- 
phonon, the electron-nonmagnetic impurity as well as the electron-electron Coulomb scattering 
into consideration. The dephasing obtained from our theory contains not only that due to the 
effective spin-flip scattering first proposed by D'yakonov and PereP [Zh. Eksp. Teor. Fiz. 60, 
1954(1971) [Sov. Phys.-JETP 38, 1053(1971)]], but also the recently proposed many-body dephasing 
due to the inhomogeneous broadening provided by the DP term [Wu, J. Supercond.:Incorp. Novel 
Mechanism 14, 245 (2001); Wu and Ning, Eur. Phys. J. B 18, 373 (2000)]. Our result agrees with 
the experiment data [Malinowski et al, Phys. Rev. B 62, 13034 (2000)] very well. We further show 
that in the case we study, the spin dephasing is dominated by the many-body effect. 

PACS numbers: 71.10.-w, 67.57.Lm, 72.25.Rb, 73.61. Ey 
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Spintronics, 1 ' 2 ' 3 where the spin degree of freedom of 
electrons is proposed to be used in place of/in addi- 
tion to the charge degree of freedom to add new fea- 
tures and functionalities in semiconductor devices, has 
attracted intensive attention recently. One of the impor- 
tant prerequisist for the realization of the device appli- 
cation is to understand the spin dephasing. Experiments 
have shown long spin lifetimes (up to 100 ns) in n-typed 
semiconductors. 4 ' 5 ' 6 ' 7 ' 8 ' 9 ' 10 ' 11 ' 12 ' 13 ' 14 ' 15 ' 16 - 17 ' 18 ' 19 In the- 
ory, three single-particle spin dephasing mechanisms, it. 
the Elliot- Yafet (EY) mechanism, 20 ' 21 the D'yakonov- 
Perel' (DP) mechanism, 22 and the Bir-Aronov-Pikus 
(BAP) mechanism 23 , have been proposed and discussed 
decades ago. 24 All these spin dephasing mechanisms re- 
quire either spin-flip (SF) or effective SF scattering. Re- 
cently, Wu proposed a many-body spin dephasing mech- 
anism where the spin conservation (SC) scattering can 
also cause spin dephasing. The new mechanism, which is 
caused by irreversibly disrupting the phases between spin 
dipoles due to the inhomogeneous broadening together 
with the SC scattering, 25,26 ' 27,28 ' 29 has long been over- 
looked. In semiconductors, the inhomogeneous broaden- 
ing can be introduced by the energy dependence of the 
g-factor 28 ' 29,30 and/or the momentum k-dependence of 
the DP term. 26 ' 27,29 In spin diffusion/ transport, it may 
also come from the different spin precession rates for dif- 
ferent moment k along the spacial gradience which we 
discussed closely in our recent papers. 31 ' 32 

It is understood that for n-typed GaAs semiconduc- 
tors, the DP mechanism is the main spin dephasing mech- 
anism at high temperature. Very recently, we performed 
a systematic investigation 33 ' 34 ' 35 of the spin dephasing 
due to the DP effect in n-typed narrow quantum wells for 
high temperatures (> 120 K) under magnetic fields in the 
Voigt configuration by constructing and numerically solv- 
ing the many-body kinetic Bloch equations. 25 ' 26 ' 27 ' 28,29 ' 36 
With all of the SC scattering (electron-electron, electron- 
phonon and electron nonmagnetic impurity scattering) 



included in our model, we are thus able to study the spin 
dephasing under various conditions. The spin dephas- 
ing time calculated from our theory contains not only 
the single particle dephasing, but also the many-body 
dephasing. Moreover, from this many-body theory, we 
are able to calculate the spin dephasing with arbitrary 
initial spin polarization (up to 100 %) which cannot be 
achieved from earlier single particle theory. Our results 
show many dramastic many-body effects for both high 
and low spin polarizations. 33 ' 34,35 Many many-body ef- 
fects we predicted, especially those under high spin polar- 
izations, have not yet been measured experimentally. In 
this communication, we will show evidence of the many- 
body effect by comparing with the existing experiments 
with small initial spin polarization. 

We start our investigation from an n-typed (100) GaAs 
quantum well (QW) with well width a. The growth di- 
rection is assumed to be z-axis. A moderate magnetic 
field B is applied along the x axis. Due to the confine- 
ment of the QW, the momentum states along z axis are 
quantized. Therefore the electron states are character- 
ized by a subband index n and a two dimensional wave 
vector k = [k x , k y ) together with a spin index a. For the 
QW with narrow width we study in the present paper, 
only the lowest subband is populated and the transition 
to the upper subbands is unimportant. With the DP 
term included, the Hamiltonian of the electrons in the 
QW takes the form: 

H = E{ £ k+[5"BB + h(k)].^) c L Ckff '+i//. (1) 

Here £k = k 2 /2m* is the energy of electron with wavevec- 
tor k and effective mass m*. a are the Pauli matrices. 
The DP term h(k) can be written as 37 ' 38 

^ IA (k)= 7 ^^-(fc z 2 )), 

h^ IA (k)= 7 ky({k 2 z )-kl), 
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hf lA (k) = 



(2) 



Here (A; 2 ) represents the average of the opera- 
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_(tt/o) 2 . 7 = 
(4/3)(m7™ CT )(l/^2m* 3 S s )(? 7 / v /l - 77/3) and 77 = 
A/ (E g + A), in which E? 9 denotes the band gap; A rep- 
resents the spin-orbit splitting of the valence band; and 
tyIqv is a constant close in magnitude to the free electron 
mass to . 39 The interaction Hamiltonian Hi is composed 
of Coulomb interaction H ee , electron-phonon interaction 
H p h, as well as electron-impurity scattering Hi. Their 
expressions can be found in textbooks. 40 ' 41 

We construct the kinetic Bloch equations by the 
nonequilibrium Green function method 40 as follows: 



Pk,<T<r' — PY,ocj' Icoh + P\,aa' Iscatt 



(3) 



Here /?k represents the single particle density matrix. 
The diagonal elements describe the electron distribu- 
tion functions pk.<x<r = /k<r- The off-diagonal elements 
Pk,i-i = Pk describe the inter-spin-band polarizations 
(coherence) of the spin coherence. 36 Note that p k ,-ii = 
p* 11 = p k . Therefore, / k +i and pk are the quantities 
to be determined from Bloch equations. 



The coherent parts of the equation of motion for the 
electron distribution function and the spin coherence arc 
given by 



at 



coh = -2<r{{gp B B + h x (k)]Imp k + /iy(k)Re/j k } + 4<rlm^ VqPk+qPk, 



(4) 



-jjf Icoh = ^9VbB + ih x (k) + hy(k)](f k i -/ k _i)+*2v q [(/ k+q i -/ k+q _i)pk-pk+q(/ k i -/k-i)]>(5) 
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where Vq = 47re 2 /[K (<z + go)] is the 2D Coulomb 
matrix element under static screening. go = 
(e 2 m*/Ko) /k=o.cr and «o is the static dielectric con- 
stant. The first term on the right hand side (RHS) 
of Eqs. (4) and (5) describes spin precession of elec- 
trons under the magnetic field B as well as the effec- 
tive magnetic field h(k) due to the DP effect. The 
rest terms are the contribution to the spin precession 
from the HF term of the electron-electron interaction. 
Pk.trcr' Iscatt stands for the scattering term. In our the- 
ory we include all of the SC scattering, i.e., the electron- 
non-magnetic impurity, electron-phonon as well as the 
electron-electron scattering. Their expressions are listed 
in detail in Refs. 33,34,35, and will not be repeated here. 
The initial conditions are taken at t = as: 

Pk|t=o = (6) 
jw|t=o = l/{cxp[(£ k - P*)/k B T} + 1} (7) 

where p a is the chemical potential for spin er. The con- 
dition pi ^ p_i gives rise to the imbalance of the elec- 
tron densities of the two spin bands. Eqs. (3) through 
(5) together with the initial conditions Eqs. (6) and (7) 
comprise the complete set of kinetic Bloch equations of 
our investigation. 

The kinetic Bloch equations form a set of nonlinear 
equations. The DP term in coherent part of the Bloch 
equation tends to drive the electron distribution out 
of equilibrium Fermi distribution, while the scattering 
term attempts to randomize electrons in k-space. There- 



fore, these two parts are coupled in the Bloch equations 
and the equations have to be solved self-consistently. 
This self-consistent solution can be performed numeri- 
cally through the method that has been laid out in de- 
tail in the previous paper on the DP mechanism in 3D 
systems. 25 By numerically solving these coupled equa- 
tions, we are thus able to study the spin procession be- 
tween spin-up and -down bands. As discussed in the 
previous papers, 36 ' 42 one can obtain the spin dephas- 
ing time as the inverse of the slope of the envelope of 
the incoherently summed spin coherence p(t) — ^ k |pk|- 
This spin dephasing time consists both the single-particle 
cffcctive-SF-scattcring-induced spin dephasing time and 
the many-body spin dephasing time. We will com- 
pare our spin dephasing time with the experiments as 
well as that predicted by the widely used single particle 
theory: 27 ' 43 ' 44 



1 _ f~dE k (f k i-f k _i)T(k) 
J™dE k {f kh -f k _ h ) ' 



(8) 



with 



T{k) = 2T 1 {k)\( 1 (kl)) 2 k*-- 1 (kl)k i 



+ 



r-\k) = / , 
Jo 



1+T 3 (fc)/Tl(fc) 

16 



7 2 fc 6 



a(E k ,6)[l-cos(n6)]d6 . 



(9) 
(10) 
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Here a(Ek,9) stands for the scattering cross-section. 

We compare the spin dcphasing times calculated from 
our many-body theory in two GaAs QW's for the well 
width a being 6 nm and 10 nm with those from the ex- 
periment by Malinowski et a/. 45 In the experiment the 
spin dephasing times for 15 and 20 nm wells have also 
been measured. However, for these wider width QW's, 
one needs more subbands than just the lowest one as 
we use in the present paper. The electron density is 



2 x 10 cm which is the background excita- 



tion measured in the experiment for both samples. The 
material parameters of GaAs for our calculation are tab- 
ulated in Table I. 46 The impurity densities Ni are taken 
to be 0.1-/V e for a = 6 nm QW and 0.05iV e for a = 10 nm 
QW. It is noted that here we use the elastic impurity 
scattering to simulate also the electron scattering by the 
interface roughness. Therefore it is understood that the 
impurity density for the 6 nm QW is larger than that for 
the 10 nm well as narrower QW suffers more scattering 
from the interface roughness. 
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TABLE I: Parameters used in the numerical calculations 

In Fig. 1 the spin dcphasing time is plotted as a func- 
tion of the temperature for 10 nm QW (a) and 6 nm QW 
(b) . The corresponding experimental data are plotted as 
dots. For comparison, the spin dcphasing time calculated 
from the single particle theory Eq. (8) from the same pa- 
rameters is also given as dashed curve in the same fig- 
ure. It is seen exactly from the figure that our results fit 
very well with the experiment data. In contrast, the spin 
dephasing time predicted by the single particle model is 
about one order of magnitude slower than the experiment 
data. Moreover, it is noticed that when the temperature 
increases, both our theory and experiment data indicate 
that the spin dephasing time slightly increases with the 
temperature whereas the single particle model predicates 
a dramatic decrease of the spin dephasing time with the 
increase of the temperature. 

These results clearly demonstrate the important role 
the many-body effect plays in the spin dcphasing. The 
single particle theory does not count for the many-body 
dcphasing due to the inhomogeneous broadening inher- 
ited in the DP term. 25 ' 26 - 27 - 28 ' 29 Therefore, it is easy to 
understand why the single particle model gives much 
slower spin dephasing time. 

Taking the inhomogeneous broadening effect into ac- 
count, one can also understand the temperature depen- 
dence of the spin dcphasing time: as the temperature 
increases, the inhomogeneous broadening is reduced as 
the electrons are distributed to the wider k-states. As 
the result, the number of electron occupation on each k 
state around the Fermi surface is reduced. Consequently 
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FIG. 1: The spin dephasing time as a function of the temper- 
ature for two GaAs quantum wells with different well width: 
(a) a=10 nm; (b) a— 6 nm. 



the inhomogeneous broadening and therefore the many- 
body spin dephasing are suppressed. It is further noted 
that this reduction is mild as a function of the temper- 
ature. Therefore, the temperature dependence is quite 
mild. 

In summary, we have shown evidence of the many-body 
effect in the spin dephasing for n-typed narrow QW's 
where the spin dephasing is dominated by the many-body 
effect. It is noted that the small well width not only lim- 
its electrons to the lowest subband, but also enhances the 
DP term as (fc 2 ) = (7r/a) 2 in Eq. (2) is greatly enhanced 
for small width and it is much larger than k 2 F for the 
density around 10 11 cm -2 . Consequently the inhomoge- 
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neous broadening is also greatly enhanced. For wide QW 
or bulk sample, the DP term is not enhanced. We notice 
from the experiment that the temperature dependence 
of the spin dephasing time is similar to what predicted 
by the single particle theory, ie., the spin dephasing time 
decreases with the temperature. A full many-body quan- 
tum kinetic calculation for the quantatative accurate spin 
dephasing time is very time consuming and is still miss- 
ing. Our earlier qualitive calculation for bulk GaAs has 



shown the decrease of the spin dephasing time with the 
temperature. 25 A thorough investigation is still underway 
to find out how large the many-body effect is for those 
samples. 
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